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The flux hne lattice of superconductors has been investigated when there exists a periodicity in 
the underlying system, such as can occur in artiflcially layered structures. For small fields parallel 
to the layers the flux lines enter the sample in sequential rows, with the possibility of jumps in 
the magnetization as new rows are created. As the held is increased these discontinuities gradually 
decrease, but there still exist transitions between states that are aligned differently to the periodic 
direction. Increasing the magnitude of the periodic potential reduces the competition between 
differently aligned lattices and tends to lock in one particular alignment. The effect of transitions 
on the shear modulus is also discussed and related to the experiments of Theunissen et al. 
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I. INTRODUCTION 

The nature of the flux line lattice in superconductors 
still attracts much interest. For an isotropic system the 
equilibrium lattice is known to be triangular. We will 
investigate situations where this may not be the case. 
In particular we shall impose a one-dimensional periodic 
structure on an infinite superconductor, and use this pe- 
riodicity to model an underlying potential. This allows 
us to model for example some of the behavior of artifi- 
cially layered structures when the field is applied parallel 
to the layers. 

The interplay between the flux line lattice and an un- 
derlying periodic potential has already been investigated 
for many different superconducting systems. These in- 
clude thin films with modulated thicknessEl, Pb-Bi alloys 
with periodic concentration of_Bi3, Pb/Ge multilayers 
with a lattice of submicron holeaj, superconducting wires 
with an hexagonal array of artificial pina3 and spatially 
modulated Josephson junctionaJ. [-. 

Using London Theory, Brongersma et aLQ modeled the 
magnetization of thin Nb-Cu multilayers with the applied 
field almost parallel to the layers. The magnetic induc- 
tion was allowed to change by permitting the size of the 
system to vary while keeping the number of flux lines 
fixed. As the applied field was increased, a series of max- 
ima in the magnetization was observed. This was due to 
the fiux lines reorganizing within the sample, where tran- 
sitions occurred between states where the flux lines form 
rows that divide the sample into equal parts. This sit- 
uation has also been investigated by numerically solving 
the time-dependent Ginzburg Landau equations coupled 
to Maxwell's equationspin a homogeneous isotropic su- 
perconducting thin filmEl. As the field applied parallel to 
the film was increased the magnetization also showed a 
series of maxima. Away from these maxima there were 
a series of discontinuities in the magnetization as more 
flux lines penetrated into the sample. 

Critical currents also show similar behavior. Disconti- 



nuities in the critical current have been observed at spe- 
cific matching fields where the numbers of rows change. 
The positions of the maxima in the critical currents is 
sensitive to the system used, with YBCO, BSCCO and 
NbCu multilayers showing different field dependencesU. 

We investigate here the properties of fiux lines within 
an infinite system, under the infiuence of a potential with 
period Ly in the y direction. We take the appfied field to 
be in the z direction. This infinite system contains a prin- 
cipal region of width Ly and approximates a layered sys- 
tem. In the absence of the potential, it is well known that 
fiux lines form a periodic lattice with hexagonal symme- 
try in an infinite system. With the periodic potential the 
fiux lines may form a lattice that is commensurate with 
the principal region and we investigate rectangular and 
centered rectangular structures with this property. The 
rectangular structure has basis vectors in the x and y di- 
rections, and the centered rectangular structure, shown 
in Fig. is characterized by the angle 2(j) between the 
basis vectors. For two different (j) the centered rectan- 
gular lattice will be the triangular lattice expected in a 
pure infinite system. These correspond to ^ = 7r/6 and 
4> = 7r/3. Fig. [|a and Fig. |l|b show cj) = 7r/6 and cj) = 7r/3 
respectively. The state with <j> = tt/G has the base of its 
equilateral triangle aligned perpendicular to the periodic 
direction. Conversely, in the state with (j) = 7r/3 the base 
is aligned parallel to the periodic direction. We distin- 
guish between these two states by referring to the cen- 
tered rectangular structure with (jj — tt/6 as the lattice 
aligned parallel to the periodic direction (the y axis), 
and the structure with (j) = tt/B as the lattice aligned 
perpendicular to the periodic direction. 

On minimizing the Gibbs free energy it is found that 
the equilibrium lattice usually has (f) only approximately 
7r/6 or 7r/3, but the above distinction regarding align- 
ments is used. The equilibrium lattice is often a competi- 
tion between these two alignments, and this competition 
is most readily seen when the periodic potential is very 
weak. 
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As the applied field is varied two distinct types of tran- 
sitions are observed between some of these different struc- 
tures. The simplest, Type A, transitions just involve the 
number of rows of flux lines within the principal region in- 
creasing by one, with no realignment of the lattice. Type 
B transitions occur between states aligned differently to 
the periodic direction, and during these transitions the 
number of rows changes significantly. 
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FIG. 1. A multilayered superconducting structure, with 
the field applied parallel to the layers. The filled circles rep- 
resent the fiux lines, which form a lattice commensurate with 
the principle region, a strip of width Ly. The lattices shown 
are centered rectangular structures, whose basis vectors are 
separated by the angle 2(j>. The figure show a) = 7r/6 and 
h) (l> = 7r/3 

The periodic potential can be chosen to favor having 
the flux lines close to cither the edges or the center of 
the principal region. This ensures one of the alignments 
of the flux line lattice is favored over the other. The 
sequence of transitions therefore can be changed by in- 
creasing the strength of the underlying potential. 

The model used to investigate the flux line lattice is 
developed in Section The formalism discussed concen- 
trates on the lowest Landau level (LLL) limit in the mean 
field approximation. Throughout we shall neglect ther- 
mal fiuctuations. The transitions are also investigated in 
the London limit. For many of these transitions notable 
jumps in the magnetization are observed. This is not the 
only property that can be effected by the periodicity of 
the system. Some of the consequences of the periodicity 
are developed in Section [V , where we investigate in par- 
ticular the shear modulus cgg. It 'oscillates' in a manner 
similar to that seen in experiments by Theunissen et al.B. 
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where = D • D and D = —ihV — 2eA, and the mag- 
netic induction B = V x A. The LLL assumes the mag- 
netic screening length is infinite. We will use the Landau 
gauge where A = Bxy, and ignore fiuctuations in the 
magnetic field. In the LLL the order parameter is ex- 
panded in the set of eigenfunctions of , using only the 
degenerate states which have the lowest eigenvalue. In 
the rectangular geometry these degenerate states are 
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where Im is the magnetic length — (<I>o/27ri?)^/^ and 
ko is the wavevector in the y-direction whose value is 
determined by the boundary conditions. The general 
LLL order parameter is 5* = Q Yl^=- 00 '^p^p^ which is 
periodic in the y-direction over a length Ly = 27r/A;o- 
Obviously states which do not have this high degree of 
periodicity of their order parameter at mean field level 
could exist for a layered system. We shall ignore them as 
even with the imposed periodicity the phase diagram at 
mean field level is extremely rich. 
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II. BASIC FORMALISM 

The experiments of Theunissen et al.i suggest the pos- 
sibility of transitions between different configurations of 
fiux lines both in the limit of low magnetic induction 
where London Theory will be valid, namely where Cge fx b 
and at higher inductions where cge oc (1 — 6^) when the 
LLL is more appropriate. We investigate the LLL first. 

The Ginzburg-Landau free energy is 



FIG. 2. The flux line lattices with rectangular and centered 
rectangular symmetry. These lattices are commensurate with 
the principal region, which is a strip of width Ly . The equilib- 
rium lattice is predominantly a centered rectangular lattice. 

In an infinite system, the equilibrium confi«|iration of 
fiux lines is known to be that of a triang lellil. In our 
problem the fiux lines are imagined to interact weakly 
with an underlying periodic potential due to the layer- 
ing. Although we ignore the direct contributions of it to 
the energy at this point, we assume that this potential 
acts to impose the periodicity Ly on the ordered state. 
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Hence, the flux lines form a lattice that is commensurate 
with the principal region and repeats over a distance ly, 
which must be a simple fraction of Ly, i.e. ly = Ly/Ny 
where Ny is an integer equal to the number of rows of flux 
lines in the principal region. We investigate rectangular 
and centered rectangular structures with this property as 
shown in Fig. |2[ The flux line lattice with rectangular 
symmetry has basis vectors li = lyy and I2 = Z^jX. The 
centered rectangular structure, also shown in Fig. ^ has 
basis vectors li = —Ix^^ + ^^y/S and I2 = Ix^ + lyfl'^, 
and is characterized by the angle 2(/) between these two 
vectors. The model used is shown for a centered rectan- 
gular structure, but is easily modified for the rectangular 
lattice. 

Within the LLL the mean-field line is defined by 
aB = where B = and as = ol{T) -\- eBfi/m. Using 
the linearized expression for a (T), a{T) — —a'(\—t) then 
as = —a' (1 — t — h) where t is the reduced temperature 
T/Tco and b is the reduced field B/Bco- T^q is the mean- 
field transition temperature while Bco is the straight line 
extrapolation of iJc2 to zero temperature. The temper- 
ature is conveniently represented by the dimensionless 
parameter ax = asiT^fiLz/ (ieBkBTY^'^ , where Lz is the 
sample thickness along the field direction. Low temper- 
atures are represented by ax — s- — 00 while high temper- 
atures correspond to ax +00. 
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FIG. 3. The dependence of Pa on the ratio R = ly/2lx. It 
should be noted that /3a{R) = I3a{1/R) 

The rescaled free energy per flux line is then 
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where Afi^x — ^o/B is the area per_flux line. The ratio 
R — ly/2lx, as used by Kleiner et al.Eil, characterizes the 
lattice. Fig. || shows the dependence of the Abrikosov 
parameter j3A on the ratio R, and it should be noted 
that fiA{R) = I3a{1/R), where (3a = m^)/^?)^). 
The square lattice corresponds to i? = 1 and the minima 
in f3A{R) occur for the equilateral triangular lattice, with 
R = 1/V3 and R = These two minima correspond 
to the two alignments of the flux lattice relative to the 
periodic direction. The state with = 7r/6 in Fig. ||a 
corresponds to the minima at i? = \/3 and the state with 
= 7r/3 has R = 1/^3. 

For a given field _B, the free energy per unit volume F^oi 
depends solely on the parameter R. If we minimize F^^i , 
the flux line lattice will be arranged in the configuration 
that minimizes Pa{R), while keeping B = 2RN^{'S?q/ Ll). 
For small fields B < ^o/Ly the equilibrium configuration 
corresponds to Ny = 1 as the smallest value of (3a{R) 
corresponds to Ny — 1. As the field is increased there 
will be transitions between configurations of different Ny 
and equal Pa- The transition between states iV^ — 1 and 
Ny = 2 occurs when B = 4<f>o/-^y- 

For larger fields the transitions occur between states 
near the minima of (3 a- Most transitions occur between 
states with i? w 1/ V3 but there are a few transitions to 
and from states near the minimum at i? = VS- AH these 
transitions occur for constant B between states of equal 
(3a{R) with a change in Ny. 



III. GIBBS FREE ENERGY 

To understand the possibilities of transitions between 
different states at constant H the relevant quantity to 
minimize is not the free energy per flux line, but rather 
the Gibbs free energy per unit volume Gvoi- The Gibbs 
free energy per flux line can be written as 
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Rewriting this in terms of the reduced fields Gibbs free 
energy per flux line becomes 
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where h ~ iiqH/Bcq. Although ax is the single inten- 
sive parameter that characterizes the free energy it is 
necessary to show the explicit dependence on the average 
induction B. Therefore, since Qvoi = QfiuxB/^oLz, at 
constant H and T, 
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where b is the re- 

duced magnetic field b = 2RNy{2'K(^ / L^) and ^ is the 
coherence length defined through the upper critical field, 
Hc^ — ^/2kHc = <i>o/27r^^. The behavior of the system 
is therefore controlled by k and 60 — 27r^^/L^. 

If the applied field h is small, the Gibbs energy 5 (/i, t, b) 
is dominated by the term (l — t— b)'^ /Pa- The minimum 
free energy flux line configuration is obtained by being 
as close as possible to having R = l/VS or R = VS. 
The range of allowed values of b is small, and as h is in- 
creased there is a series of first order transitions at which 
the lattice reconstructs. As the lattice undergoes a re- 
construction between different states 'if{Ny,^/3R) there 
are jumps in the magnetic induction b. Except for the 
transition ^'(1, 1) ^(1, 3) this involves a change in the 
number of 'rows' of fiux lines and these reconstructions 
can be classified into two categories. The simplest involve 
just an increase in the number of rows for two states with 
R « l/-\/3 or R fa V3- Some of these Type A transi- 
tions include ^-(2,1) ^ "^(3,1) and 4'(4, 1) ^ ^'(5,1) 
and occur between states with the same alignment to 
the periodic direction. The other set of transitions in- 
volve a change in the alignment of the flux line lattice. 
The Type B transitions occur between states near the 
different minima of /3a (-R) and involve either a reduction 
in the number of rows combined with a large increase in 
R, or an increase in the number of rows and a reduction 
in R. These transitions include ^'(3,1) ^'(2, 3) and 
^1^(4, 3) — > ^(7,1). By putting in ascending order the 
values of KNy where K = 1/ \/3 or a/3 we can see the 
appropriate sequence of states. This sequence of states 
looks similar to that in the previous section, but the 
states possess a smaller range of values of R. Transi- 
tions occur between states with different (3a, different b 
and different Ny. 
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FIG. 4. A representation of how R changes as H increases. 
The range of values of R that minimize the Gibbs free energy 
are shown within the shaded region. 

For very large applied fields, such that 2k^(& — h)^ is 
the dominant term in g(h,t^b), there is a wide range of 
equilibrium values of R that are obtained as h is varied. 
The lattice distorts from a triangular lattice to ensure 
b K h. However, transitions between states with different 
Ny still occur, but these transitions involve very small 
changes in the magnetic induction b. As the applied field 



h is increased there is a smooth crossover between these 
two limits as the two terms in g{h, t, b) compete with each 
other. This results in the allowed values of R spreading 
out from the two values of 1/ VS and a/S- 

These transitions can be shown schematically by look- 
ing at the allowed values of R that minimize the Gibbs 
free energy. The transitions corresponding to small fields 
h are between states defined by R being very close to 
either l/\/3 or V^S. As h increases, the range of values of 
R that minimize the free energy increases, and this leads 
to the fractional jump in b dropping. 
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FIG. 5. The jumps in the magnetization vary as the ap- 
plied field is increased, decaying to zero in the limit h ^ 1; 
for Ly = lOA, hi — 50. The points indicate the change in the 
magnetization when the potential is very weak. The heavy 
line is a guide to the eye, whose physical significance is dis- 
cussed in the next section. 

The size of the jumps in the magnetic induction asso- 
ciated with these transitions is given in Fig. |. These 
jumps increase at small h as approximately ft. 2 but de- 
cay to zero in the limit h 1 following the heavy curve 
in Fig. ^. The inset describes in detail the initial tran- 
sitions, showing the points lying on or below the start 
of the heavy line. Within the LLL approximation the 
fiux lines are not excluded from the sample at small h as 
i/ci = 0. The flux lines always form a centered rectangu- 
lar structure and we only consider behavior for magnetic 
inductions larger than i?mm- This is the magnetic in- 
duction that corresponds to a simple triangular lattice 
with Ny = 1 and R = l/VS. The first point in the in- 
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set corresponds to the transition between states ^'(1, 1) 
and ^(1,3). The subsequent sequence of these centered 
rectangular structures is ^'(2, 1), ^'(3, 1), ^'(2,3), ^'(4, 1), 
'I' (5,1) and is easily predicted to the limit h ^ 1. The 
Type A transitions, such as 5'(2, 1) v[/(3^ i) (Point 
3) and ^'(4,1) *(5,1) (Point 6), lie very close to 
the heavy line. All the other points in the inset de- 
scribe Type B transitions, *(1, 1) ^1^(1, 3) (Point 1), 
*(1,3) -> *(2,1) (Point 2), 4'(3, 1) ^ ^-(2,3) (Point 4), 
^-(2, 3) ^-(4, 1) (Point 5), and ^-(5, 1) ^ *(3, 3) (Point 
7). The Type B transitions tend to lie well below the 
heavy line, and these transitions involve a large change 
in the number of rows. The heavy line describes the sit- 
uation when Type A transitions only occur, which can 
happen when the strength of the potential is increased. 
This is discussed in Section 

The transitions can also be observed as the tempera- 
ture is varied. For small iVj, these 'melting' transitions 
only occur over a very small range of the applied field h. 
As the range of allowed values of R increases the range 
of field over which the lattice can be seen to melt also 
increases. Fig. ^ shows the transition between states 
^'(8,> 1) and ^-(9, < 1). 
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A. London Theory 

The LLL will not be strictly applicable in the limit of 
small b. The LLL approximation may be valid only in the 
region where dAb/dh < in Fig. |. By increasing bg it is 
possible to shift the maximum in Ab{h) to smaller values 
of Ny but London Theory should be used to investigated 
the properties of the system with small values of Ny. 

The London free energy is 



^Lon ^l^JZ jj dr?l^a/3(r. - r, )dr; 



(7) 



where Vap (r) is the potential defining the local magnetic 
induction B{r)^ — / drfVa/sir — r^), the parameters 
{a, (3} correspond to the {x, y, z} components, and {z, j} 
sum over the contributions from the all flux lines. We 
investigate an isotropic system, where the Fourier trans- 
form of the London potential for straight flux lines is 
VapCk) — Sai3S{'k)/{l + A^/c^), A being the London pen- 
etration depth. The cutoff function S'(k) = exp(— ^-^fc^) 
removes the divergences within London Theory due to 
the absence of the flux line cores. The use of differ- 
ent cutoffs—within London Theory has been discussed 
previousljitj. 

We use similar geometries to the previous section and 
minimize Gloji = ^Lon — B ■ H with respect to the po- 
sitions of the flux lines. Again we assume the flux lines 
form a lattice commensurate with the principal region. 

The low field behavior using the LLL was dominated 
by the minimum in Pa^R)- It might be thought that 
using London Theory could remove the multiplicity of 
transitions. 

The equilibrium lattice is always a centered rectangu- 
lar structure, and the flux lines do not initially enter the 
principal region in a single straight row. This is different 
from a thin slab geometry, where straight rows of flux 
lines (rectangular symmetry) are frequeii^^ the equilib- 
rium structure due to the surface barrier 
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FIG. 6. The 'melting' transition between states with 
Ny = 8 and Ny = 9 ioi h = 6.9 x lO""*, Bo = 2 x 10"^ 
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FIG. 7. The change in the magnetization as a function of 
the applied field h for Ly = lOA and k = 50. The upper heavy 
line shows the jumps associated with the LLL. 



V{y) = Vo(cosh(a?/) — 1) where —Ly/2 < y < Ly/2 and 
assume the potential is periodic as described above. The 
potential is zero in the center of the principal region, but 
depending on the sign of Vq will be either attract or re- 
pel flux lines form the edges of the principal region. This 
potential then favors states with R = 1/ ^/S or those with 
R = -s/S. Assuming the potential is still weak, we include 
just the first order correction to the ground 

state energy. 

The potential favors the Type A transitions between 
states near the same minimum of (3a{R) and the inclu- 
sion of this potential reduces the number of Type B tran- 
sitions. As Vq is increased it is then possible to remove 
transitions between the two minima in I3a{R) and for 
transitions to occur between states with R ^ I/a/S or 
between states with R ~ -s/S- In the figure analogous to 
Fig. 1^ the allowed values of R would be one shaded band 
or the other and there would be no transitions between 
them. 

The repulsive (attractive) potential ensures that the 
lattice chooses alignments aligned parallel (perpendicu- 
lar) to the periodic direction, i.e. the y-direction. The 
jumps in the magnetic induction follow the heavy curve 
in Fig. ^. This behavior also occurs in the London limit, 
where the inclusion of a potential again favors the ex- 
istence of states permanently aligned either parallel or 
perpendicular to the y-axis. 



The flux lines enter as a centered rectangular struc- 
ture with Ny — 1, and there is a smooth crossover from 
^(1,< 1) to *(!,> 3). The first transition occurs be- 
tween states ^'(1, > 3) and \l/(2, < 1), and the sequence 
of structures observed is then the same as the LLL limit. 
At each of the transitions there is an associated jump 
in the magnetic induction, but Fig. ^ shows that these 
are smaller than in the LLL limit. Again, during the 
Type B transitions (point 3, ^'(3, 1) ^'(2,3) and point 
4, ^'(2,3) — > ^'(3,1) ) the change in the magnetiza- 
tion is much smaller than during the Type A transitions 
(point 2, ^-(2, 1) ^ ^-(3, 1) and point 5 *(4, 1) ^'(5, 1) 
nearby Therefore, despite the absence of the pronounced 
minimum in the Gibbs free energy the flux line lattice 
shows similar behavior to the LLL limit and the same 
sequence of states is observed. 



IV. NON-ZERO POTENTIAL 

In the previous section, the potential has not only 
been assumed to be periodic, i.e. V(y) = V{y + Ly) 
but also to only restrict the choice of suitable flux line 
lattices. This may be either an extremely weak poten- 
tial or one whose contribution to the total energy is 
zero. Now we investigate the equilibrium lattice when 
the periodic potential does contribute to the overall en- 
ergy. We define the potential only in the principal region 




FIG. 8. The changes in the magnetization for Ly 
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showing both the LLL limit (filled squares) and the London 
limit (open circles). 
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The behavior in Fig. ^ for Ly/\ = 10 is controlled 
by Bq. Reducing the width of the periodic potential at 
fixed penetration depth A, the maximum in A6(/i) moves 
slightly, to smaller Ny. Fig. || and Fig. || show the jumps 
in the magnetization for different widths. The LLL (filled 
squares) and the London limit (open circles) are shown 
together and can be compared with Fig. || and Fig. |^. 
The heavy lines indicate the transitions that occur in the 
strong potential limit, and all transitions disappear in 
the limit h ^ 1. The behavior is very similar for the 
three different widths shown, but increasing the width 
increases the number of transitions observed. 
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FIG. 9. The changes in the magnetization for Ly = vTOA 
showing both the LLL limit (filled squares) and the London 
limit (open circles). 

Theunissen et ali have investigated the properties of a 
NbN/NbsGe bilayer under the influence of an applied 
current and magnetic field. A regular array of chan- 
nels was etched through the thin NbN top layer into the 
NbaGe. The NbN layer, with its much higher critical cur- 
rent than the NbsGe layer effectively acted as a pinning 
center for the flux lattice. By measuring the shearing 
force needed for the flux lines to move along the chan- 
nels, the shear modulus of the flux line lattice could be 
measured. The shear modulus shows a characteristic be- 
havior C66 fx. 6(1 — 6^) as the field is increased. In the Lon- 
don limit C66 c>c b while in the LLL limit cge c>c (1 — 6^). 
Superimposed on this functional form was an oscillatory 
function that reflected the finite nature of the channel 
width. 



The shear modulus cge oc d'^Fyoi/da^ where a is the 
shear angle. In the LLL limit, Fyoi oc The 
Abrikosov parameter can be written as a sum over all 
reciprocal lattice vectors Q, Pa — X]q ^xp — Q^/2C/^, 
where U is the inverse magnetic length U"^ = 2'K/Aunit, 
A-unit being the area of the unit cell. This allows direct 
calculation of the shear modulus for all lattices. 
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FIG. 10. The shear modulus for the commensurate flux line 
lattice with no added potential, showing the oscillatory nature 
as a function of the magnetic induction, for Ly — -\/TOA and 
K = 50. 

The modulus changes dramatically as the lattice re- 
constructs. Fig. |l^ shows Cge in the LLL , where the os- 
cillations characterize the rearrangement of the flux line 
lattice within the channel. This shows more oscillations 
than the experiments of Theunissen et alu, and the intro- 
duction of the (attractive) potential models the experi- 
ments more closely. The size|-and positions of the tran- 
sitions that Theunissen et al.cl observed indicate transi- 
tions only with AiVy = -1-1 and for R ~ -\/3- Therefore, 
the effects of the interaction of the flux lines with the 
underlying periodic potential are likely to be important, 
although the interaction with impurities may hinder ex- 
tensive reconstructions with large changes in Ny. 

V. CONCLUSION 

For an inflnite system it is well known the flux lines 
form a periodic lattice with hexagonal symmetry. Once 
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there is competition with an underlying potential this is 
not necessarily the situation. The flux lines still favors 
lattices close to this ideal lattice, and transitions can be 
seen between different structures. Using both London 
Theory and the LLL approximation we have investigated 
some of the properties associated with these structures. 
As the applied field is increased there are notable jumps 
in the magnetic induction, and hence in the magnetiza- 
tion and critical current. These transitions occur in two 
ways. The Type A transitions just involve an increase 
in the number of rows of flux lines in the principal re- 
gion. Type B transitions occur between states aligned 
parallel and perpendicular to the periodic direction. The 
different alignments of the lattice correspond to the two 
competing minima in Pa{R)- Increasing the strength of 
the periodic potential ensures one of these alignments 
is favored over the other and eventually only the Type 
A transitions increasing the number of rows occur. For 
Vq less(greater) than zero, the potential favors the flux 
lines being at the edges (center) and the flux line lattice is 
aligned parallel(perpendicular) to the periodic direction. 

The transitions are most easily seen as the applied field 
is increased. However they also occur as the temperature 
is varied and the lattice can appear to melt as the tem- 
perature is increased. Also associated with these tran- 
sitions are changes in the physical properties of the flux 
line lattice, such as the shear modulus. The matching 
fields and the nature, of the cge indicate the experiments 
of Theunissen et aLcl are best modeled by an attractive 
potential, which pulls the flux lines to the edges of the 
region. 
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